Skip to main content
Internet Archive's 25th Anniversary Logo

Full text of "Microlensing in phase space II: Correlations analysis"

See other formats

Mon. Not. R. Astron. Soc. OOO.H1IT7IC2006') Printed 5 February 2008 (MN I*TeX style file v2.2) 

Microlensing in phase space II: Correlations analysis 

A.V. Tuntsov* & G.F. Lewis f 

A29, School of Physics, University of Sydney, NSW 2006, Australia 











Accepted 2006 April 6. Received 2006 April 4; in original form 2006 January 16 


Applications of the phase space approach to the calculation of the microlensing auto- 
correlation function are presented. The continuous propagation equation for a random 
star field with a Gaussian velocity distribution is solved in the leading non-trivial 
approximation using the perturbation technique. It is shown that microlensing modu- 
lations can be important in the interpretation of optical and shorter-wavelength light 
curves of pulsars, power spectra of active galactic nuclei and coherence estimates for 
quasi-periodic oscillations of dwarf novae and low-mass X-ray binaries. Extra scatter 
in the brightness of type la supernovae due to gravitational microlensing is shown to 
be of order up to 0.2™ depending on the extent of the light curves. 

Key words: gravitational lensing - methods : statistical 
stars : oscillations - galaxies : active - dark matter 

methods : analytical 


Statistical approaches to gravitational microlensing have 
a long his tory th at extends back to at lea st the mid 
1980s (e.g.JPaczvnski 1986, Kavser. Refsdal fc Stabell IQSl 

IPeeuchi fc Watsonlll987fl when it was realized that in gen- 
eral, a large number of stars contribute to the deflection 
field, and as a result, it is nearly impossible to model in- 
dividual micro lensing light curves in detail unambiguously 
JSchneider fc Weiss 1987: Katz. Balbus fc Paczvhski 1986) . 
Since then, a number of methods have been developed to cal- 
culate the statistical properties of microlensing light curves. 
Apart from a few early p apers, the methods employed are 
numerical - ray-tracing ([K avser, Refsdal fc StabcU 198d 


Wambseanss. Witt fc Schneider. .1992 : Schneider fc Weisd 
1983) and image-tracking techniques JLewis et alJ Il993 


Witdll99a) - because they alone are capable of calculat- 

ing the magnification probability distribution function over 
a large range of magnification values. 

However, it was soon realized that analytic methods 
could help solve the complementary problem of correla- 
tions in the magnification values at nearby positions in 
a microlensing ligh t curv e. Based in part on the work of 
[Peguchi fc WatsonI lll987l) . ISeitz fc Schneided lll994) deveb 
oped a method to calculate the deflection angle probability 
distribution function and successfully applied it in a deter- 
mination of the flux autocorrelation function for a source 
observed through a single lens plane for a variety of the pa- 
rameters iSeitz. Wambseanss fc SchneideJll994l) . However, 
this method did not attract much further interest because of 

* E-mail: tyomichOphysics. 
t E-mail: 

substantial technical difficulties in its implementation. Nein- 
dorf (2003) took a further major step forward by developing 
a simple approximation to the correlated deflection proba- 
bility distribution that could be applied to many interesting 
astrophysical situations. 

In a companion paper (Tuntsov fc Lewis 2006, hereafter 
Paper I), we presented a similar approach in which a phase 
space (or intensity) description of the deflection process is 
used before projecting intensities onto observable quantities. 
The method makes clear some of the important approxima- 
tions used in studies of gravitational lensing and is easily 
generalized to include many lens planes along the line of 
sight, which is especially interesting in the case of a cosmo- 
logically distributed population of lenses. To avoid arbitrary 
'slicing' of the space into individual lens planes, we have also 
proposed a continuous approximation to the propagation of 
intensity moments arguing that for average quantities there 
is no conflict between the thin lens approximation and a 3D 
distribution of lenses; we call this 'nearly' 3D lensing. 

In Paper I we have also solved the problem of the prop- 
agation of flrst-order moments. However, this is of little in- 
terest by itself as the results follow trivially from the conser- 
vation laws and the required symmetries of the initial condi- 
tions. In the present paper we apply these results to the non- 
trivial problem of the propagation of second-order moments 
of the flux; traditionally, the most interesting projections are 
the autocorrelation function and its Fourier transform - the 
power spectrum. 

In section 2, we review the results of Paper I in relation 
to this subject and present a general model for the initial 
conditions. In section 3, the deflection potential is calculated 
in the flrst non-trivial approximation for a random star field 
with a Gaussian distribution of velocities. The formal so- 

Tuntsov & Lewis 

lution of the propagation equation using the perturbation 
technique is obtained in section 4, which is followed by the 
calculation of the flux power spectrum for a point- like source 
in the leading order of the deflector density. For clarity, most 
of the calculations of this section are presented in the Ap- 

Section 5 discusses the application of these results. First 
we consider modifications to the power spectra of periodic 
and quasi-periodic sources. We show that known microlens- 
ing populations of the Galaxy can introduce ~ O.f^-deep 
modulations in the pulse profiles of a kpc-distant pulsar and 
significantly affect coherence estimates for quasi-periodic os- 
cillations of X-ray binaries and dwarf novae in the mHz to 
Hz frequency range. Then, corrections to the autocorrelation 
functions of aperiodic sources are considered; in particular, 
we find that in the 0.01 — 0.1 range of redshifts used for the 
calibration of type la supernovae brightness, the extra scat- 
ter due to gravitational lensing can amount to 0.04™ — 0.2™ 
depending on the light curve sampling time-scales. We con- 
clude with a short list of refinements to the method that can 
increase the quality of our predictions. 


In paper I we have considered the transformation of the 
statistical moments of the intensity field as light propa- 
gates through a universe filled with gravitational lenses. The 
second-order moments are defined as the average value of 
the product of intensities at a pair of positions in the phase 
space 7V/(^ 9 (ri, r2, ai, 0:2, ti,i2), defined by the spatial po- 
sition of a light ray in the lensing plane r^, its direction a.i 
and the moment of observation ti, {i = 1,2): 

M2 (ri,r2,ai,Q;2,ti,i2) = E J (ri, ai,ii) J (r2, 02, ^2) • (1) 

The expectation is taken with respect to an ensemble of 
realization of the intensity field; further discussion of this 
ensemble is given below. 

It was shown in Paper I, that the conservation of surface 
brightness in gravitational lensing - that is, the conserva- 
tion of the intensity along the ray - results in the linearity 
of transformation on M2 as the light rays propagate from 
one lens plane to another and get defiected at them. It is 
therefore convenient to consider M2 as a vector in some lin- 
ear function space V (M^\. In particular, M2 can be repre- 
sented in any other functional basis; certain properties of the 
transformations often suggest particularly convenient bases. 

Under the thin lens approximation it is possible to in- 
troduce two operators, termed deflection "I? and focusing 
F, which describe how the intensity moments are modified 
when the light rays are broken on the lens planes ($) and 
when they propagate freely between the planes {F). Then, 
we derived a Fokker-Planck-type equation to describe the 
continuous limit of this treatment, and to answer the ques- 
tion of how the intensity moments behave as a function of 
depth D along the line of sight: 





Here d<i>/d_D and d-F/d_D are effectively the generators of <!> 
and F, describing the changes in M2 due to these operators 
applied over an infinitesimally small distance dD. 

In the case of the second-order moments, the focusing 
operator in the above equation remains the same: 

dF _ _ >r^ 
df5 ^^ 2^ 




+ 2^a, 


Unlike the first-order case, there is a difference in the 
symmetry properties of the deflection operator. The thin 
lens and small angles approximations still apply, and there- 
fore the deflection probability depends only on the deflection 
angle /3i = a'^ — cti, and not on a;; this leaves the associated 
deflection operator diagonal (proportional to a (5-function) 
in Ti coordinates - Fourier conjugates to cti. In the focusing 
operator, we simply replace cti with —id-n- 

At the same time, deflection probabilities are no longer 
entirely independent of temporal or spatial coordinates as 
the deflection angles and potential time delays of the two 
rays are correlated. Dependence on t means that the asso- 
ciated defiection operator is no longer diagonal in the tem- 
poral frequency basis. However, recalling that the potential 
time delays turned to be unimportant in the observationally 
accessible range of frequencies, we will ignore them in the 
subsequent analysis and in doing so avoid a difficult integral 
term in the equation for M2 . Then, the defiection probability 
operator is diagonal in ri,Ti,ti basis. 

In the stationary and spatially uniform case, probability 
densities for the light rays to be deflected by (3i depend on 
the difference between the spatial and temporal coordinates 
associated with the deflections: 

<l(ri,Ti,ti,r2,r2,i2 -> ri, ri,ii,r'2, T2, ^3) = (4) 

'I>(ri,r2,ti,i2,Ti,r2)(5(Ti - ti)S(t2 - T2) 
x5{r', - ri)5(r'2 - r2)3{t[ - U)5(t'2 - ^2), 

so the equation for the second moment assumes the following 

dDM2 (D;ri,r2,ii,t2,Ti,T2) = 

/ 1 ' 1 


1 + TT^ii^t; 

2^. .1^.1 ' 2c "^*=^ 

+P2 (ri,r2,ti,i2,Ti,T2) JM2 (_D;ri,r2,ti,t2,Ti,r2) 

where p2 is the derivative of $ with respect to D. 

To complete the statement of the problem, one also 
needs to specify the initial conditions M2|d=o = ^2- In 
the present study, we will restrict ourselves to the case con- 
sidered in Paper I: an isotropic source located at D = —Ds- 
Then, for the first moment one has: 

Mi'(k,a;,r)=/(k,a;)— -exp 

'"^ (r + kD,)^ 





where f{\i.,u)) is the Fourier transform of the source profile. 
However, there is an important difference between the 
first and second order moments. In Paper I the ensemble 
over which one should average to obtain the moments was 
not fully defined because it served as a conceptual tool in the 
general derivation. Two alternatives were provided for the 
definition - one either looks at independent patches of the 
sky and averages quantities of interest along different lines 
of sight or observes a certain source long enough to consider 

Microlensing in phase space II: Correlations analysis 3 

different patches of the Ught curve as independent elements 
of the ensemble. 

In this paper, we are more concerned with practical ap- 
plications and for this reason we will stick to the latter defi- 
nition as it is closer to the situation encountered in observa- 
tions. In this case, an internally consistent definition of the 
averages can only be provided for sources that are statisti- 
cally stationary. In the first order, this leads to the trivial 
result - in the absence of a special point in time, all modes 
with (jj 7^ have zero amplitude and stay so because of the 
homogeneity of equations, while the u — mode is preserved 
due to conservation of energy; even for the alternative defi- 
nition using different sky patches, it is hard to think of any 
reason for a certain time position to stand out, as noticed 
in Paper I. 

This can be expressed in a formal way. Indeed, shifting 
the origin of the time axis by a constant T results in a phase 
factor e'""^ in front of the Fourier transform of a function 

iujT , 




the same applies to all averages. Therefore, in a statistically 
stationary case, when EtF — EF for all T one has 

ErFiuj) = e'""EF(a;) =EF{iu), 


Hence, EF{uj) = for all tj / 0. 

Similarly, for the second order moments one has 

ETF2iLOi,u;2) = e'^'^'+'^^^^EF2{uji,iU2), 



and the only non-trivial statistically stationary components 
are those with a;2 ~ — i^i = i^/2; we presume w ^ from 
now on. When F{t) itself is a real quantity - as is the case 
for intensity and flux - F{—lj) = F*{uj), and the function 

Pf{uj) = EF(-lj/2)F{lu/2) = E \F{Lu/2)f (10) 

is the power spectrum of F oscillations and represents a 
Fourier transform of the autocorrelation function of F. 

If one now assumes that the source brightens and dims 
synchronously over its entire spatial profile so that it can be 
separated into a temporal and spatial parts: 

/(k,<.^) = <7(k)/iH 


the initial conditions for the second moment of intensity 
distribution of an isotropic stationary source located at _D = 
—Ds can be written as: 




X exp 


((n + kiD, 

5(ki)ff(k2) (12) 


where Q — {uji+uj2)/2, lo — lu2—uji. The average g{'ki)g{'k2) 
is used because correlations may be present in the source 
profile. From now on we will drop the factor 5(0.) which 
arises from the infinite extent of time axis and thus can be 
considered as a constant through regularization. 


In calculating the deflection potential p2 we will again 
assume that the distribution of the deflectors is statisti- 
cally uniform and stationary and ignore any correlations 

between individual deflectors; we will also assume Poisso- 
nian distribution for the number of deflectors in a lensing 
plane. In t his case, the Markov-Holtsmark-Chandrasekhar 
iChandrasekharnl943i : jNeindort, 20031 method gives the fol- 
lowing result: 

P2 (Ti,T2,r,t) = iy\\C\\ (p(Ti,T2,r,t) - 1) , 


where f is the spatial density of the deflectors, r = r2 — ri , 
t = t2 — ti, \\£\\ is the area of the lensing plane, and p is the 
Fourier transform of the defiection probability density with 
respect to the two deflection angles, calculated for a single 
deflector in the plane. 

The probability density itself can be easily obtained 
since, for a single deflector of mass m, the deflection an- 
gles uniquely determine its positions pi,2 with respect to 
the rays: pi — Vi— mfii/ 01. Therefore, for a uniform plane, 


p(/3i,/32,r,i) = 



(j) \ m 



Here the fraction is the Jacobian of the pi,2 -^ /3i,2 trans- 
form, while <j!>(p, t) is the probability density for the position 
change p over a time interval t - essentially, the velocity 
distribution of the deflectors. The simplest choice will be 
to assume a Gaussian distribution of these velocities with 
average v and dispersion g: 

r /_ ..-in2" 

(p,t) = Vl-Ko t j exp 




However, Fourier transforming 11411 even in this sim- 
plest case is a very hard task. We can, however, expand it in 
Taylor series in ti,2 and restrict ourselves to the most im- 
portant contribution in this expansion. This is a necessary 
exercise given that the problem for l|^ is in ten dimensions - 
numerical methods will not be easy to implement here unless 
we have some analytic understanding of the solutions. 

From the normalization of the probability density, the 
zeroth-order term in Taylor expansion of p(ti , T2 , r, t) is 
unity and it cancels out in 11311 . The coefficients in front of 
Ti and T2 are the average defiections of the two rays which 
equal zero for a uniform lensing plane. Therefore, the first 
term of the expansion is of the second order: 



{{{t, ■ f3^f) + {{t2 ■ 02)')) (16) 


+ o{ti,T2] 

where the angled brackets correspond to averaging over the 
distribution 11411 . 

For symmetry reasons, the averages on the flrst line do 
not depend on r and t, and therefore they cannot depend 
on the direction of ri,2 either. Thus, 

(((n . /3i)^) + {{t2 ■ (32)'}) = i{/3^) (rf + r^) . (17) 

The other average in H16^ is responsible for the correlations. 
It can be rewritten in a simpler form 

2{(ti ■/3i)(t2 •/32))(r,i) = {P'}t^Co ir,t)T2, (18) 

where Co is the correlation matrix 

C;Mr,t) = 2{/3^)-^{/3i/3^)(r,i) (19) 

The average in the angled brackets can be found using the 
joint distribution function of /3i,2: 

Tuntsov & Lewis 





Later we will need the the Fourier transform of the correla- 
tion matrix as a function of r and t: 

2 ,, exp [i (k ■ r + cjt)] 





For a Gaussian distribution of velocities HSU this integral 
can be done analytically. One first integrates with respect 
to r (or equivalently, r — ?nAl//3) over an infinite plane 
to obtain a Gaussian-type integral with respect to t, while 
the integrals with respect to /3i,2 decouple. As a result, one 
obtains a degenerate matrix: 

where a marks the direction of k: k = kea and 

2 m^ 

C(k,w) = ^/---^exp 






where we have assumed that the upper limit m/po on (3 is 
determined by the lens physical size po while the lower one, 
po/R comes from the extent of the lensing plane R — Apo- 
For relativistic compact objects, such as black holes and 
neutron stars, po/m ~ 0(1) while for white dwarfs po/m ~ 

In the case of continuous distribution of the deflecting 
surface density E(p, f), the deflection angle power spectrum 
can be related directly to the statistical properties of E. For 
a spatially and temporally uniform case one has (E is in 
mass per length squared units): 


4G^ fc'fe^ 



C* fc4 

where Pe is the projected density power spectrum, which 
can be found using an appropriate form of the Limber's 
equation. In the formula above we have neglected retarda- 
tion effects, therefore it is valid when the motion of deflecting 
matter is slow compared to the speed of light. 

However, the smooth-matter description requires care- 
ful treatment of the correlations along the line of sight and 
will not be considered in this paper. 


Neglecting the higher-order terms in Hlfcil , we may write JSJ 




+ ( 19x1 ■ 9ri + 19x2 ■ 9r2 + TT^^l ^«1 + TT^l^ ^t^ ] ^2 




In the absence of correlations represented by the first line of 
the right-hand side a dense system of eigenvectors Ma (and 
corresponding eigenvalues \a) of the remaining operator is 
provided by various products and associated sums of the 

M4ri,r2,Ti,T2,ti,i2) = MW(ri,Ti,ti)Mi^'(r2,T2,i2) (26) 

Aa — y^ai + ^a 



with a composite index a = (01,02), where MaJ , XaJ are 
the solutions for the first-order equation: 


X'a'/M^l> (r, T, t) = [id^ ■ dr + y/^dtj M^J' (r, r, t) (28) 


As was found in Paper I, these represent harmonic modes in 
r and t multiplied by an angular functions Cmni''') depending 
on the spatial and temporal frequencies k, uj as parameters 
(see Appendix IXt. 

The functions Ma^ form a complete orthonormal ba- 
sis of the solutions, while their products of form 12611 do 
the same job for the problem considered in this paper. It 
is therefore convenient to rewrite 12511 in the basis of Ma- 
Namely, let 

M2(D;ri,2,Ti,2,ii,2) = Vcae^°-°M4ri,2,Ti,2,ii,2) (29) 

with a — (ki,2,<jJi,2, ?wi,2,'T.i,2) running through all possible 
values. When there are no correlations a set of Ca — const 
chosen such as to represent the initial conditions at D = 
gives a solution of 12511 . 

In the presence of correlations, which we will treat as a 
perturbation, we allow these 'constant coefficients' to vary 
with D to obtain the following system of linear equations in 
the perturbed case: 

Ca — / ^aa'^ 

(A ,-\a)D 


where the dot denotes differentiation with respect to D. The 
matrix elements Caa' are defined as 

Caa' = -M\{(3'} {MlT^CoT2Ma') 


The perturbation method then proceeds by introducing 
a dummy parameter e to characterize the order of perturba- 
tion eC and expanding Ca{D) in its power series: 

Ca{D) = c°(Z?) + €cl{D) + e\l{D) + . 


such that (^{D) functions correspond to no correlations at 
all, c\ (D) are proportional to the first power of perturbation 
and so on. A sequence of successive corrections to Ca{D) can 
now be obtained by equating coefficients at common powers 
of e: 

Ca = 0, 


(x^,-\a)D n 



and so on. 

From the first of the above equations cJJ = const and if 
we choose these constants to represent the initial conditions 
M2 at D = 0: c° = (m],M2), the initial conditions in aU 
other orders will be uniform 

4(0) = 0, 



Physically, the zeroth-order term describes just the partial 
loss of the phase coherence in different modes as a result of 
geometric time delays due to uncorrelated random deffec- 
tions, as explained in Paper I. 

Microlensing in phase space II: Correlations analysis 5 

Further, we restrict ourselves to the first non-trivial 
term in the expansion 13211 . It is proportional to the strength 
of correlations in the deflection angle at different points and 
describes how these correlations translate across to the in- 
tensity field. By mixing the intensity modes of one- ray statis- 
tics CaiCaJ exp[(Aai -|- Aaa)-^] in a correlated way, extra cor- 
relations in the combined field are created, which can be 
described by solving 1341 using c° = const: 

e c„(D) = > c„, 

dD' Ca'e 

XaiD~D') X iD' 


The power in the first-order term is also turned into extra 
signal by correlated deflections, but the higher-order contri- 
butions will be neglected in this study. 

Although the intensity field provides a complete de- 
scription in the sense discussed above, the intensity is rarely 
observed in the real life due to the limited angular resolution 
of our telescopes. The quantity that is normally measured 
is the fiux F. 

_F is a linear functional on the intensity field and, just 
like intensity, defined for every spatial and time coordinate 
r, i. Thus, in calculating it, one only needs to consider the 
projection in the angular sector - namely, consider t — 0. 
We can therefore prescind from the exact parametrization of 
r and t dependence - e.g., the temporal Fourier harmonics 
of the flux observed at points ri,2 

-F2(ri,r2, 1.^1,(^2) = (37) 

will be those of the intensity projected at t = without any 
additional transformation. 

In the approximation we use, there are two contribu- 
tions to these quantities, associated with the zeroth and 
first-order terms in the perturbation theory expansion. They 
both depend linearly on the initial conditions and there- 
fore one can introduce two D-dependent functionals Fj ' £ 
V''{M^) that calculate the flux Fourier components as func- 
tions of M2 at respective accuracy. 

In the Appendix]^ these functionals are calculated with 
respect to Ma functions. Since they form an eigenbasis of the 
unperturbed problem, F° is diagonal in this representation. 
This is not the case for F^: correlated deflections can pro- 
duce the power in the flux Fourier spectrum at frequencies 
where the flux power spectrum due to initial conditions is 

A stationary, isotropic and synchronous source is de- 
fined by just two functions: the average source spatial profile 
Fourier transform g(ki)g(k2) and fluctuations power spec- 
trum P{u)). Therefore the expectation value of the autocor- 
relation function Fourier transform (ki,2 = K =F k/2): 

F°'^ (K, k, u^;D)=Y. i^:,}(23)P(a;')<?(k'i)<;(k'2) 


with d — (K,k,aj); here we use the summation instead of 
the integration sign to avoid unnecessarily long notation. 

The coefficients F^^, (D) are given in the Appendix ^ 
Eis an explicit expression for F^^, and as an integral of an 
explicit expression for F^^,. They describe the fraction of 
variability power that, in the presence of correlated defiec- 

tions on path length of D, is transferred from the original 
(a;', K', k')-mode to a mode at (cj,K,k). As noted above, 
F^ is diagonal in this basis. 

However, a question that is more interesting from the 
observational point of view concerns the value of the fiux 
temporal variability power JT observed at a certain point R 
(or autocorrelation between a pair of points ri_2 = R=F i'/2) 
in the plane of the observer. This can be obtained by sum- 
mation of all amplitudes multiplied by their Fourier mode 
value at the given point (s). 

^"•^(^,0)= / dLj'd'kid'k2d^k'id k2F^'^'(D) 



,, exp[-i(ki • ri) - i(k2 • r2)] 


For the non-trivial term, coefficients F^^, depend not only 
on ij,ui' but spatial wave numbers as well, and therefore 
the summation can only be done for a certain source profile 

For the purpose of illustrating this method we choose 
the simplest case of a point-like source at the origin of the 
source plane and a single observer at the origin of the plane 
of the observer ri,2 = 0. To facilitate comparison with the 
results of Paper I, we assume that the source is located at 
D — —Ds ^ and that D > space is filled with deflectors 
of constant number density u. We perform the calculations 
in Appendix |n| here only the necessary results are present. 

In order to visualize the results we consider a 'line pro- 
file broadening' function ^{lj,uj' ,D) that represents broad- 
ening of a single line in the power spectrum at a certain u': 
P{cl>) = P{iu')5{iL!—uj'). 4> is defined as the ratio between 1391 
and the value of the flux power spectrum T^{u',D)\ 
which would have been observed at D were the lensing pop- 
ulation absent. The latter can be obtained by noticing that 
the power spectrum is quadratic in flux which obeys the 
inverse- squares law when the light rays are not deflected. 
For the above choice of the source profile and the position 
of the observer (see the Appendix!^ 

.F°(^',g)^ ./^""l. ^0(v^^ 

{Ds + DY 


and therefore 

i^°'^{Lo,uj',D) = {Ds + Df /d'Kd^K'd'kd'k'i?°^J(D)(41) 

The part of JT" due to uncorrelated deflections is of the 
second order in the deflectors' density v. As shown in Pa- 
per I, in all relevant astrophysical situations the density of 
deflectors is very low in terms of a dimensionless quantity 
rnD^JvLo/c that controls the strength of gravitational lens- 
ing by a three-dimensional distribution of compact lenses. 

The dominant non-trivial contribution to ^j) in this 
regime is supplied by the term due to correlated deflections 
J-^ ~ Olyy/v). Assuming that parameters u,m,a,v are con- 
stant along the line of sight one can introduce the frequency 
difference scale 




6 Tuntsov & Lewis 

an amplitude parameter 




m 8n\nAv {D^+Dy 

and write down the following approximation for i/) (we drop 
superscript '1' because this is a dominant non-trivial contri- 
bution) : 

12 exp |"-«V2cr2] 
7/;(Ac^, D) ^ ^ i * ' 


/27rlnAAa;o KAuJo'Ds+D'a 

where a scaled version of the line profile function 

-,g (44) 


q I dxw {s,xq) exp 




\(J X J 

An auxiliary function it)(s,^) is defined in the Appendix IHl 
and proportional to ^~^ for ^ 2> 1. Since A 2> 1 we can safely 
neglect the second Bessel function in the square brackets 
when evaluating (1451 except for some extremely low values 
of y when it does induce a small correction; Jo can also be 
neglected in subsequent formulae 15UII and 1521 . 

Figure displays the behaviour of the line profile ^ as 
a function of the frequency difference variable y for s = 
and w/o" = 1. It has a flat central part and then falls off 
gradually. The behaviour of '^ is rather hard to study ana- 
lytically as it has an essential singularity in the integrand; 
nevertheless, it can be understood qualitatively by noting 
that exp[— y^/2a;^] effectively cuts out x < y from the in- 
tegration region and using known form of w(^) and Bessel 

There are two different regimes for <&(?/). At s = 0, for 
q < 0(1) the height of the flat part is proportional to q^ 
and it extends to y ~ q~^ followed by a '^{y) ~ ci~^y~^ tail 


q V 




For q > 0(1) the height of the central part is oc g and 
the extent is again 0{q~^); however, the fall-off rate here is 
up to J/ ~ C(l) where it turns to the same 


<? ^y~ 


■^ q y 
^ tafl: 


z y t 

0{q-^)<y<0{l) , g<0(l).(47) 

y . 

Oscillations of the Bessel functions enter the region not dom- 
inated by the exponential fall-off at different values of y 
which results in a few 'wiggles' in the resulting 'i!{y) curve 
seen in Figure 

In the opposite case of a distant source viewed through 
a localized overdensity of lenses, 1 — s ^ 1, the behaviour is 
generally the same but q parameter should be replaced with 
its scaled version 

q =q{l-- s 


m / Stt In A vD^ 


Figure 1. The line profile function ^{y, s,v/a,q) as a function 
of frequency argument y for s = 0, v/cr = 1 and a number of 
different values of the amplitude parameter q. The profile has a 
fiat central part of the height oc q (for q > 0{1)) or q^ (otherwise) 

extending to yc ^ 0{q 

up to 

y ~ C(l) and as q~^y~^ afterwards. Certain 'wiggles' can be seen 
in the shape of ^ which stem from the behaviour of the Bessel 
function in 1451 . 

the scaling arises from the fact that the main contribution 
to w(s,^) comes, for ^ ^ 1, from the region of i such that 



Different cuts through 'i/{y,q) can be obtained by con- 
sidering ^I" as a function of q, as shown in Figure |5| here, 
again, we assume that the source is inside the deflectors- 
filled volume: s = and choose v/a = 1. For a given fre- 
quency difference parameter y, <!/ first increases {y is a part 
of the central plateau) as either <& ex g^ for q < 0(1) or 
^ oc g for g > 0(1), and then decreases as q~^ . The enve- 
lope curve of the ^(g) family thus represents the height of 
the central part as a function of the amplitude parameter 
g. Taking into account the scaling (1481 . s 7^ (to be more 
precise, 1 — s ^ 1) curves can be obtained by a shift of the 
ones given in Figure |5| to the right. 

For completeness, let us also state the dependence of 
the line profile on the velocity parameter v/a. At y <ti 1 
^ is nearly independent of v/a while at the shoulder and 
tail it scales as roughly exp[v'^ /2a^]. As a result, the shoul- 
der and tail of the product exp[—v^/2a^]'^{y,s,v/a,q) are 
insensitive to the velocity while the centre of the profile is 
significantly suppressed for v/a 3> 1 and gently connected 
to the shoulder by a practically featureless curve. 

The total fraction of the initial power added to the spec- 
trum by correlated deflections can be obtained by integra- 
tion of tp with respect to Alu: 

y{D)= / dAujiP{ALu,D) 


12exp[-«V4a2] ^ f v^ \ ,, f 
i^A ^°(,4^j^l 


Ds + D' 



Y{s, q) = Jd^ ^w {s, [Jo i^/q) - Jo (AC/g)]' (50) 

Microlensing in phase space II: Correlations analysis 7 

- - Y(q) 

., Y,(q) 


,,, Y^(q)/qY(q) 


0.01 0.1 1 10 100 1000 10000 

Figure 2. The line profile function ^{y, s,v/a,q) as a function 
of the amplitude parameter q for s = 0, v/a = 1 and a number 
of different values of y. For every frequency argument value y, 
function ^(q) reaches a maximum at qdy)'- Vciqc) = J/ and then 
falls off as y enters the wing of the line profile. The envelope 
curve of this family represents the value of * near the edge of the 
central part ydq)- 

Figure 3. Functions Y{s,q), Yi{s,q), Yi{s,q)/q and the dimen- 
sionlcss effective 'line width' Yi (s, q)/qY{s, q) as a function of the 
amplitude parameter q for the source embedded in the deflectors- 
filled volume: s = 0. Y and Yi both saturate for large values 
g 3> 1 and have oc q tails for g <§; 1. The dimensionless effec- 
tive 'line width', on the contrary, is a monotonically decreasing 
function of q. 

Another quantity of interest is the effective 'line width' 
defined as: 

(Auj) (D) = y'^ dAiu \Auj\ i/. {Aiu, D) 



^ /T exp [^V4a^] y,(s,g) 

'^"V2 hiv^/^o^) qY{s,q) 

Y,{s,q) = / dUw{s,0 [Jo (C/g) - Jo (AC/g)] = 


we need to make use of the absolute value of Alu in such a 
definition because as can be seen from the asymptotics of 
w{s, ^) and Bessel function the second moment of Aui does 
not exist. 

Using II45I I one finds that, for q < 0(1), the total frac- 
tion of the power created by correlated deflections y is pro- 
portional to q and is divided in roughly equal parts between 
the central plateau and y~^ tail. They also make approxi- 
mately equal contributions to the integral with respect to y 
in 1511 1. which is independent of q. As a result, the 'effective 
line width' (Ao;) oc q~^ . 

For q > 0(1), we find, using 1471 . that y saturates for 
g 3> 1; most of the power is split roughly equally between the 
central plateau that is high (oc q) and narrow (oc q~^), and 
the y~^ shoulder that extends out to Alo ~ 0(A(x;o), only 
a small fraction ~ 0{q~^) is in the y~'^ tail. The 'effective 
line width' in this case (Aw) oc In q/q. 

Figure |3 shows the results of numerical integration 
in (I50II52II . They agree well with the approximations above. 

An important property of the broadening function i/; in 
the low-density approximation is that it only depends on 
the frequency difference Au = u' — uo. The convolution the- 
orem then makes the transformation of the autocorrelation 
function (which, by virtue of the Wiener-Khinchine theorem 
is the Fourier transform of the power spectrum) very sim- 

ple: one just multiplies the autocorrelation function by the 
Fourier transform of ^ and adds the result to the original 
autocorrelation function. 


When discussing the application of the results we have ob- 
tained, it is important to interpret the results of observations 
correctly. Within the model we have used, all real observa- 
tional programs have two fundamental limitations. 

Firstly, it is the finite temporal resolution of observa- 
tions T, making the information on variability at high fre- 
quencies inaccessible to observations. For the power spectra, 
this can be modelled by multiplying them by some window 
function of the frequency which decays rapidly for u) > t~^ 
and does not affect it for low frequencies. A good choice here 
could be a Gaussian normalized such that its value at zero 
is unity: 

Wt{lli) = exp 


Another important limitation is the finite extent T of any 
real observational program. This brings about an inability 
to estimate variations on very long scales and results in a 
power spectrum where adjacent frequencies are no longer 
independent; as a result, fine spectral features are smeared 
out. This can be modelled by convolving the power spec- 
trum with a similar window function of width ~ T~^; since 
the total power should stay constant, the function must be 
normalized such that its integral is unity. An example of 
such a function is an integral-normalized Gaussian 

WTii^,l-^ ) 






Let us now turn to some specific examples of the ob- 
jects and study how their variability is affected by gravita- 

8 Tuntsov & Lewis 

tional microlensing. We will first consider periodic or nearly 
periodic sources, sketching the kind of effect gravitational 
microlensing brings about for objects observed through the 
microlensing population of the Galaxy on distances of order 
1 kpc; it will be convenient to use the frequency domain in 
this study. 

Then, aperiodic sources will be considered with a spe- 
cific application to the extra scatter in the flux of stan- 
dard candles; here the language of autocorrelation functions 
will be more appropriate. An important point regarding the 
standard candles is that they are most interesting on the 
cosmological scales while the analysis present in this work 
is essentially Newtonian in its use of space and time. In ad- 
dition, the parameters of the lensing population on these 
scales can only be guessed; therefore the results obtained in 
cosmological applications should be viewed cautiously but 
do give a general feeling of what the effect of microlensing 
on these scales should be. 

In both cases optical (and shorter) wavelengths will be 
assumed, because at radio wavele ngths interstellar scintilla- 
tions are much more effective (e.g.. lScalo fc Elmegreerl2004i 
FLazio et al. 20041. 

5.1 Pulsars 

The most prominent features in the power spectra are dis- 
played by objects which are periodic or close to being pe- 
riodic. One of the most interesting for us will be pulsars. 
The way their variability is affected is an interesting ques- 
tion on its own; in addition, they will serve as a tool to help 
us interpret the action of spectral broadening on this simple 

The power spectra of periodic objects such as pulsars 
display an array of spectral lines at multiples of the spin fre- 
quency ojp; the relative heights of these lines depend on the 
pulse shape. Broadening of every line is a typical feature of 
the multiplicative correction with a slowly- varying function. 

Indeed, let us consider a periodic source of flux 

J^p(i) = ^yl„e'""''*, A„ = yl! 


with a power spectrum 

Ppiuj) = 2^ \^n\ S{u) ~ ncjp), 



modulated by some function H-e(t), with e varying so slowly 
that its Fourier transform e(Aaj) is essentially contained in 
the region Auj <^ uip. 

Then, for the power spectrum of Fp{t) — Fp{t){l + e(i)) 
one has^ 

Ppiio) = Pp{Lo) + Y^ \A„f P,{lo - nujp 


+ 2_^ A„A^e{uj — nu}p)e*(i 


LU — mUJr, 


where Pe is the power spectrum of e: Pe(Aaj) = |e(ALi;)|'^. 
The last line in l|57|l can be neglected due to assumption 

^ We drop the infinitesimal correction to the discrete part of the 

above, and comparing what is left to the results of our cal- 
culation in the previous section, we can interpret ^(Ao;) as 
the power spectrum of the magnification factor. 

Let us estimate the quantities involved for a typical 
pulsar at a distance of 1 kpc observed through the pop- 
ulation of white dwarfs of mass 0.5 Mq with the space 
density of 0.01 pc~'^ (as recent studies along different lines 
suggest - HolbcrE, Oswah & Sion '2002"; Hafizi ot al. '200^; 
.Qppenheimer et al. 2001 : Han fc Gould 2003') : assuming 
they have halo kinematics, we may choose o ~ 250km/s 
(e.g., iBinnev fc Trema ind Il98 q). For such w hite dwarfs, 
on/m » 3 X 10^ ijShapiro fc Texikolskvlll98,'Jl . Given that 
pulsars generally have high spatial velocities, a choice of 
v/a — 2 seems to be reasonable. 

Then, for InA = 30, the amplitude parameter becomes 
g « 17 and the frequency scale Aluq = 0.05 s~^. Given a 
relatively high value of the velocity parameter, tp starts off 
as an increasing function of Ao; as the latter changes from 
zero to Aluc ~ 0.01 s~^, where it turns to a shoulder followed 
by the y~^ tail for Alo > 0.3 s~^; a small step is present at 

The total extra power added to a line by gravitational 
microlensing is, according to 14911 . y « 15% which is con- 
tributed mostly by the shoulder - the central part accounts 
for about a tenth of this power while the tail makes up a 
fraction of order few percent. The tail contributes slightly 
more to the integral 1521 1. and the effective line width, as 
given by ||MJ, is (Acu) ^ 0.06 s~\ 

These frequencies scale according to 1G14II with param- 
eter a - e.g., the frequency difference scale would be twice as 
large if the defiectors' velocity dispersion was only a half of 
the one assumed and twice as low if the defiector population 
was two times faster. 

Our predictions for the effect of gravitational microlens- 
ing can be interpreted in a simple way for pulsars with spin 
frequencies greater than about ls~^. Namely, we suggest 
that the pulse height is modulated by a slowly (and, pre- 
sumably, randomly) varying function with an rms of 0.08 
stellar magnitudes on the time-scales T > Auj~^ ~ 100 s; 
the time resolution required r < 1 s. 

The actual pulsars show a degree of pulse-to-pulse vari- 
ability that is much stronger than the expected 0.08™ con- 
tribution from lensing. However, unlike the lensing signal, 
these variations seem to be uncorrelated even on interpulse 
scale and therefore can be suppressed by averaging a large 
number of pulses. This does not necessarily mean a long 
time-scale - the lensing signal correlation time is of order 
seconds. Therefore, ceteris paribus, the lensing signal is more 
likely to be detected in high-frequency pulsars. 

5.2 Quasi-periodic and related oscillations 

Another phenomenon of considerable interest in relation 
to the subject of this paper is quasi-periodic oscillations 
(QPOs) generally defined as broad - up to of order the cen- 
tral frequency - lines seen in the power spec tra of various 
astrophysical sources (e.g.. Ivan der Klial2005^ : more coher- 
ent features, with higher centroid frequency to width ratios 
Q, are known by a variety of names, deflned phenomenolog- 
ically and often specific to a class of sources or even to the 
particular state of a given source. 

Within the Galaxy, these oscillations are associated 

Microlensing in phase space II: Correlations analysis 9 

with accretion onto compact objects - black holes, neu- 
tron stars and white dwarfs - from low-mass normal com- 
ponents. The frequencies of these oscillations span a wide 
range from mHz for cataclysmic variable stars that con- 
tain accreting white dwarfs to Hz and even kHz for low- 
mass X-ray binaries where the accretion proceeds onto neu- 
tron stars or black holes. Normally, features with Q in the 
range of up to a hundred or so are called QPOs, while 
for oscillations with greater coherence values terms such as 
normal, horizontal or flare branch oscill ations (N/H/FBOl 
in the case of low-mass X-ray binaries (Ivan der Klial2005i 
iHasinger fc van der Klij[l989l) and dwarf nova oscillations 
(DNOs), including longer-period ones (IpDNOs), in the case 
of cataclysm i c variables are in use IIKluzniak et alj l2005i 
IWarnejEooi : IWarner fc WoudlJl2004) . For convenience, we 
use the short term QPO more generally to mean all of the 
above when necessary. 

It is believed that careful modelling of these spectral 
features can hold potential clues to the physical processes 
in the inner part of the accretion discs and that QPOs may 
become important probes of the strong gravitational and/or 
magnetic flelds of the centra l compact objects (Psaltis 2001; 
iHeUieil Eooi : iMauchd l20o3) . Therefore, one has to inter- 
pret these features properly, including accounting for pos- 
sible contamination of the signal en route to the observer; 
obviously, gravitational microlensing, in the form described 
above, is one such possible 'contaminant' able to signifl- 
cantly modify the intrinsic power spectrum - although it 
does not change the the centroid frequencies of QPOs, the 
total power in the features and their width can be consider- 
ably affected by lensing. 

When making an estimate of the effect of lensing on 
the power spectra of these sources, the only change that 
seems necessary for the parameters of the lensing problem 
given in the previous subsection concerns the velocity - a 
choice of v/a = 1 is probably more sensible in this case. 
As a result, the central part is no longer suppressed but 
otherwise tp looks nearly the same. The total extra power 
added by lensing is slightly greater at y ^ 0.25, and the 
effective line width is lower: (Alo) ~ 0.035 s~^. 

We can therefore conclude that gravitational microlens- 
ing has the potential to influence the total power observed in 
a given feature of the power spectrum and should be taken 
into account in their interpretation. With regards to coher- 
ence estimates, lensing makes a negligible contribution to 
the line width of kHz QPOs but cannot be ignored in the 
interpretation of spectral lines at lower frequencies and may 
be the dominant effect in broadening lines at mHz frequen- 

These results can be summarized in a graph that 
loosely classifies oscillations of different kinds in frequency- 
coherence coordinates (Figure [IJ. Apparently, whenever the 
lensing effective line width is of order 2-kv/Q or greater and 
y ^ 0{1), it should be taken into account. Therefore, lens- 
ing becomes important for Q > 27ri^/{Atj) - the regions that 
lie above the straight lines in the plot. 

5.3 Aperiodic sources and standard candles 

Sources that do not show any prominent features in their 
power spectra, are more conveniently discussed in terms of 
the autocorrelation function. As mentioned in the previous 

Figure 4. A very rough classification of different spectral lines 
in terms of the central frequency u and coherence parameter 
Q = u/Au (Ai/ is the spectral line width). The left part of the 
plot is occupied by different features seen in the power spectra of 
cataclysmic variables while on the right one can sec more rapid os- 
cillations typical for accretion onto relativistic components. Grav- 
itational lensing becomes important in the interpretation of co- 
herence when broadening line width {Aui) becomes comparable to 
u/Q - these regions lie above the straight lines seen in the figure 
which are plotted for g = 5 X 10^, q = 17 and q = 0.5 correspond- 
ing to the distances of 10 kpc, 1 kpc and 100 pc, respectively; the 
other parameters are fixed as in the text. The total extra power 
corresponding to these three values of q are 0.32, 0.25 and 0.01 - 
that is, broadening contributes a little in the latter case although 
the line is formally quite broad. 

section, tf}{uj,uj') only depends on the frequency difference 
Alu = ij' — u! in the approximation we use. As a result, the 
correction Af{t) to the autocorrelation function f{t) can be 
found by multiplying f{t) by the Fourier transform of ^: 

A/(t) = /(t)V'W- 


However, we should take into account the fact men- 
tioned in the beginning of this section - namely, limited 
temporal resolution of observations and thus the inability 
to observe high frequency signal. To model this, we multi- 
ply ■tl){A(jj) by the function Wt{i^) defined by (I53I I and then 
Fourier-transform the result to obtain an 'observationally- 
corrected' version of ip (we include \/27r factor in the defini- 
tion of ip): 

■tPr{t) = / dije-'"^'Wr{co)i){uj). 


Both tp and Wt are even functions oi uj, so the exponent in 
the last expression can be replaced with the cosine function. 
Denoting to = Aug^ we can rewrite the above formula 
in a form, similar to 14411 : 

Mt) - i^^-PJ-f/^-^l^, / * D 

In A 

to' Ds+D' a' 


*r (?7,s,w,g) 






This integral also has to be done numerically. We have 
found the following rearrangement of the integral to pro- 

10 Tuntsov & Lewis 






^v ^ 

^v \ 





S. \ 


^- 10- 



\ \ 







10 ^ 





Figure 5. The autocorrelation correction ^^ as a function of 
rj = t/to for q = 10, s = 0, tj/cr = 1 and a few values of the time 
resolution parameter r/qto — zero, one, ten and a hundred. It has 
a flat central part followed by ri~^ tail. 

Figure 6. The height of the central plateau ^T\ri=o as a function 
of the time resolution parameter r/to for q equal to one, 0.1, 0.01 
and 0.001; the velocity parameter ti/cr = 1 and s = 0. The quan- 
tity is weakly sensitive to the time resolution for t < (Aa;)~^ ^ 
toq/lnq and decreases oct^^ after that. 

duce fast and stable convergence for various values of the 
parameters (again, we assume A is so large that the Bessel 
function it appears in can be neglected): 



dx X wi^s, x) Jq (x/q) j{rjx/q, rAujQx/q, v) 

with an auxiliary function 


dy cos hy exp 

1 + i 2 

— ?; — y 

loivy). (63) 

At u = 0, q is simply a Gaussian of width \/T+t? while at 
non-zero v, this function displays a few oscillations around 
zero value before being exponentially suppressed on the 
same scale ~ O (\/l + t^j . 

Based on the similarity between j{h) and a Gaussian, 
one can expect that ^t has an effective width of order 
(Aa;)~^, which is approximately ~ toq/lnq for g > 1. Simi- 
larly, we expect that r should assume the role of this char- 
acteristic time when it exceeds the latter; cutting off some 
power at higher frequencies also leads to certain suppression 
in the value of ^t at low t. 

Figure confirms these expectations. It shows the be- 
haviour of ^I/t as a function of r; = t/to for q — 10 and a 
few values of the parameter r/io - 0, 10, 10^ and 10"^ that 
correspond to t{Alo) of zero, approximately one, ten and 
a hundred. One can see a flat central plateau extending to 
approximately the greater of (Aa;)~^ and r, followed by a 
ri~^ tail which is nearly independent of r. The height of the 
plateau as a function of r/to for a number of q values is 
shown in Figure 1^1 

In Figure |51 a small value of the velocity parameter 
v/a = 1 is assumed. For higher velocities the Fourier integral 
for $ is dominated by non-zero frequency modes, and ^t 
drops below zero outside the flat central part, then tends to 

zero value in the same, rj 
Figure |7| illustrates this. 

', manner from the negative side. 




tAm /q=3 

,,, TA(O|j/q=10 




















Figure 7. Behaviour of the autocorrelation correction <&,- as a 
function of r] for a high value of the velocity parameter v/a = 5 
and a number of different values of the time resolution parameter 
r/qto — one, three and ten; the amplitude parameter g = 10 and 
s = 0. To incorporate both the different scales and the negatively- 
valued region, the linear scale is used but the functions are divided 
by their values at r; = 0. 

Although not easy, it is interesting to compare these 
findings to the results of previous works on calculation of 
the lensing au t ocorre lation function. Both Seitz et al. (1994) 
and iNeindor J (|200.'l ) perform these calculations for a two- 
dimensional distribution of lenses, and the resulting curve 
depends on a number of parameters, including the effective 
convergence and shear and scaled radius of the source, while 
this work assumes a (nearly) three-dimensional distribution 
of lenses with no external shear, and calculations are per- 
formed for a point-like source. 

However, inspection of Figures 9 and 11 of Seitz et al. 
(1994) shows that their autocorrelation function r{Y) does 
drop off as approximately an inverse square of Y in the 
shearless and parallel-to-shear cases for Y > 3 while 

Microlensing in phase space 11: Correlations analysis 11 

perpendicular-to-shear curve is clearly not a power law; it 
should be noted, however, that this impression is based on 
measurements of the actually published plots. On the con- 
trar y, the growth of the autocorrelation function in Figure 
1 of lNeindora i2003r seems to be shallower than the inverse 
squares; however, it is the result of a best-fit procedure to 
the actual autocorrelation function observed in one of the 

Let us now describe how an autocorrelation function of 
a distant source is affected by the transformation we have 
presented. For any source, the intrinsic flux autocorrelation 
function fit) should tend to the square of the flux average 
value for large i. Let us for simplicity assume this value to 
be unity: fit) —* 1 for t — ^ cxd. On top of that, one has all 
the information about second-order statistics of variability, 
normally presenting itself as a decreasing - not necessar- 
ily monotonically - function of t. In particular, the t — 
value of the autocorrelation function gives the flux disper- 
sion: cr^ = /(O) - 1. 

When / is subjected to the transformation described in 
this paper, the correction is proportional to /(i), not the 
variable part /(t) = f{t) — 1, and as a result, even the con- 
stant brightness sources are affected: 

fit) = fit) - 1 = m + m (1 + m) ■ (64) 

As can be seen from 16011 and the results of the calcula- 
tions in this section (e.g.. Figure |HJ, t/j is less than unity for 
any realistic In A. Therefore, for mildly variable sources with 
/(O) = cr'^ < 1, the above formula can be rewritten in the 
most simple form 



- that is, the microlensing correction simply adds up to the 
intrinsic autocorrelation function. 

The amount added is actually rather modest. For In A = 
30, the factor in front of exp[— u^/2o-'^]^ in 1601 is approx- 
imately 0.41 and ^(0) saturates at about 1.3 for z <^ 1. 
As a result, ijj{t) < 0.5, and as with the intrinsic fit), the 
maximum is at t = 0. 

For sources with some memory of their past, one expects 
fit) to be smooth at f = and therefore the '(/'(t) plateau 
is hard to spot. Whether gravitational lensing becomes a 
dominant contributor to /'(i) therefore depends largely on 
the shape of the intrinsic function fit). As ip cc t~^ at large 
t, it necessarily becomes dominant if fit) decreases as a 
power law f~" with n > 2 or more rapidly - e.g., exponen- 
tially. The well-sampled observed autocorrelation functions 
are indeed often fit with exponentials but t his is mostly for 
the purpose of simp licity and convenience llHawkinj|l996i 
Ide Vries et alj|2004f) . A simple transformation of the auto- 
correlation function - the structure function - is more often 
used to characterize variability, but despite its simplicity, it 
is non-linear on /(f) and therefore is less suitable for our 

The power spectra of various quasars and AGNs have 
also been studied o bservationally in a number of works 
fGiveon et al." '1999'; 'P altani et all Il997l : iFiore et alj Il998l: 
[Edelson fc Kirpal 19951) . Traditionally, their variability is 
divided into a short and long time-scale; the former has 
stronger amplitudes and is thought to arise from the pro- 
cesses in the inner engines of quasars while the origin of 
the long-term variability is less clear and a variety of differ- 

ent mechanisms have been proposed (e.g.. lBregmanlll99d : 
IZackrisson et al.ll2003r : it is quite likely, that a number of 
mechanisms actually contribute and microlensing could be 
just one of those. We therefore are only concerned with the 
long-term variability. 

According to the theory presented in this study, any 
strong signature in the power spectrum of a quasar will be 
broadened by a universal function i/)(Aaj). In particular, if 
a quasar or AGN has only one dominating line at tj = 0, 
its entire power spectrum will be given by this function. 
According to estimates given in the previous section, for 
q ^ I most of the power resides in the Auj~^ shoulder of 
tp. The observed power spectra on long time-scales are well 
described by a featureless power law with and index ~ —2 
outside the central region where the spectrum flattens (it 
is hard to estimate it here, though, due to the limited ex- 
tent of the light curves) . Edelson & Kiroal ( 1999|) give an 
estimate of —1.74 ± 0.12 for the slope of the X-ray fluctua- 
tions power spectrum of a nearby Seyfert galaxy NGC 3516; 
the characteristic frequency of the break in the spectru m is 
ujc « 27r/l month ^ 3 x 10"*^ s~\ Iciveon et alJ lll999ll re- 
port an average slope of —2.0 ± 0.6 (at frequencies range 
10"^ s~'^ <uj < 10"'' s"^) for the sample of ~ 40 Palomar- 
Green quasars observed regularly in the optical over seven 
years; based on the simulated sensitivity of their method, 
they conclude that they cannot reliably infer these quanti- 
ties for individual quasars in the sample. They do, however, 
report two cases with an apparent break in the power-law 
slope at a frequency cuc ~ Sx 10~^ s^^. However, estimation 
of the power law index is still difficult with the present-day 
data accuracy and coverage. The total variability power in 
these spectra are of order a few to about twenty percent 

Therefore, although it would be extremely preliminary 
to claim that these data support our predictions, they 
certainly do not rule out the microlensing hypothesis at 
the present time. Some other aspects of this variability, 
such as chromaticity and (tentative) asymmetry of light 
curves, however, arc more difficult to accomodate within it 
(Zackrisson et al. 2003). Nevertheless, it would be interest- 
ing to consider whether the microlensing-induced broaden- 
ing, if ever detected, can place meaningful constraints on 
the global microlensing population. The model presented in 
the previous section is probably too crude to describe the 
actual distribution of microlensing bodies in a cosmological 
situation with strong clustering of lenses, and more work 
needs to be done to infer quantitative characteristics of this 
population from the available data. 

Even if microlensing never shows up as a dominant com- 
ponent of the autocorrelation function or the power spec- 
trum of a source, it still contributes significantly to the dis- 
persion of the flux. Indeed, *I'(0) saturates very quickly with 
g > 1 (see Figures 13 and |SJ and for the lensing population 
considered in this section, this happens relatively quickly, 
at distances of order a few kpc. As a result, for moderate 
values of v/a, one gets a top-up to the magnitude of order 
0.1™ — 0.2™, as estimated in subsection 5.1. 

In the language of frequencies, it is clear that this vari- 
ability is spread over a range of timescales from roughly to 
to qto and consequently can be detected only if the time 
resolution r of observations is better than to; a fraction of 
order r/iqto) of this variability will be 'smoothed out' if 

12 Tuntsov & Lewis 

r > to. The larger time-scale T is responsible for the sample 
completeness - for shorter light curves large variations are 
likely to be left outside its coverage; this is not an issue if 
the dispersion of an ensemble of similar sources is consid- 
ered, which is the case when we calculate the extra scatter 
in the brightness of standard candles due to microlensing is. 

Perhaps, the most interesting of such questions deals 
currently with the magnitude scatter of distant type la su- 
pernovae used in cosmology to measure the luminosity dis- 
tance as a function of redshift for constraining the evo- 
lution of the Universe JRiess et alJIlOQa : IPerlmutter et alJ 
Il999r . The full lensing aspect of th e problem is beyond 
the scope of this paper (see, e.g.. iHolz fc L indcr 2005; 
lAmanuUah. Mortsell fc GoobaiJIJOol iHui fc Green..2005., for 
recent studies of the role lensing plays here) as the geometry 
of the Universe assumed here is flat and static and cannot 
describe the high-redshift Universe. 

However, the treatment of this paper still applies to the 
Universe at moderate redshifts up to at least 0.01 - 0.1, and 
this is precisely the redshift range used to determine the 
zero-point of the rcdshift-magnitude relation and its uncer- 
tainty ( Wood-Vascy ct al. 2004) ; an estimate of the effect 
in cosmological settings should also give some idea of what 
kind of time-scales are involved. 

The cosmological microlensing population is still quite 
a mystery. For the economy of the hypotheses we will still 
assume that the deflectors are unseen white dwarfs of 0.5Mq 
which are uniformly distributed with density v. It can be 
written in terms of the fraction of the critical density in 
deflectors fid and Hubble constant Ho as: 

_ SHo^d 
which leads to 






where the Hubble distance Dh = c/Ho has been introduced; 
one can also recognize yrnD as the radius of the Einstein- 
Ch wolson ring for a lens of mass m seen from a distance 
D llSchneider. Ehlers fc Faicolll992l) . D/Dh is very close to 
the redshift z for small values of the latter in a flat Uni- 
verse. For a c urrently preferred value Dh 
JLahav fc Lid dlc 2004), one has 

^"'"^ %o/10'^km/s 






An estimate fid In A ~ 1 can be regarded conserva- 
tive since we know that MACHOs contribute around ten to 
twenty per cent of t he dark matter in ou r neighbourhood 
jAlcock et alJ |200C| iLasserre et alJ l2000l : ISadouled Il999t 
lEvans fc Belokurov^l2004^ ■ and the dark m atter share in the 
comp osition of the Universe is Q,dm ~ 0.27 JLahav fc Liddld 
12004^ . As a result, g ^ 1 at redshifts z ~ 10~^ and above - 
which includes all intergalactic scales. It also seems sensible 
to choose the same velocity dispersion a = 250km/s (e.g., 
iTormen et al]ll993l: iBahcaU fc Siandll996ll and the velocity 
parameter v/a = 1. 

For these parameters, V'(O) _„ ~ y ~ 0.35. The time- 
scales are to = ^'-^n = 20 s and 

gio ~ 8 X 10*^ s 

po/lO'J cm 






At z = 0.01 the latter time-scale is approximately 10 days. 
The brightness of a supernova is normally derived from a 
fit to the obs erved light curve with the extent of order ten 
to fifty days JRiess et alJl2004r . which, therefore, should be 
considered the time resolution. Given that Vt(0) oc r~^ at 
these scales, we can expect the extra scatter in the magni- 
tude to be 


2 -0.4 In 10 




at these redshifts. Such an extra scatter is comparable to 
the intrinsic scatter in the brightness of type la supernovae, 
which is currently estim ated at am ~ 0.15"* iJTonrv et alJ 
1200,4 IAs 

et al.ll200Rll . and therefore should be taken into 

At greater z ~ 0.05, which is roughly the target redshif t 
of the SuperNova Factory project ijWood-Vasev et al.l2004^ ■ 
the larger timescale becomes toq ~ 10^ days and the whole 
extent of the high-quality light curve typically used in esti- 
mating its brightness t < 50 days is not sufficient to signifl- 
cally alter the scatter due to gravitational lensing. Therefore, 
this project should be able to put some constraints on the 
spatial density of lenses v. In these settings, we cannot ex- 
pect these constraints to be particularly useful because most 
of crjji"'^ comes from the first ~ 10 — 100 kpc of the optical 
path. If a method could be devised to actually estimate the 
autocorrelation function from these data, it has the potential 
to serve such a probe. However, as mentioned above, other 
more stable objects seen to high redshifts (e.g., quasars) 
could be better targets for such studies. The numerical pre- 
dictions for the amplitude of the autocorrelation function 
correction and the time-scales involved are, of course, the 
same as for the supernovae as long as both can be consid- 
ered point-like. 

5.4 Conclusions 

In the present paper we have calculated the autocorrelation 
function of the flux of a point-like source through a pop- 
ulation of uniformly distributed point-like gravitational mi- 
crolenses using the lowest nontrivial order approximation for 
the deflection potential defined by 11311 . In the last section 
we have also identified several potential applications in illus- 
tration of the method. Although, given the simplicity of the 
approximation, the estimates we obtained should be taken 
rather cautiously, they do show that the effect is important 
and cannot be ignored. 

It is sensible to conclude by recapping briefly the major 
assumptions used in this study, which should give us the idea 
of its applicability and show how and in which directions it 
can be improved. The 'physical' approximations employed 
here are the use of the geometric optics and thin lens ap- 
proximation as well as neglecting the potential time delays 
when writing down the propagation equation 0. The first 
of these was shown at the very beginning to be an excel- 
lent one in most situations in the gravitational lensing (e.g., 
Schneider et al. 1992), the other two were extensively dis- 
cussed in Paper I and Section 2 of the present paper. 

The major mathematical approximations are the trun- 

Microlensing in phase space 11: Correlations analysis 13 

cation of Taylor series in 11311 on the first non-trivial, corre- 
lation, term and the perturbation technique used in the solu- 
tion of the propagation equation. The justification for series 
truncation comes from the fact that coefficients in front of 
the higher-order powers are all proportional to the appro- 
priate powers of the characteristic deflection angles, which 
only become of order unity very close to the lens; this never 
occurs for the deflectors considered in this paper; however, if 
one wants to consider black holes as deflectors, it will be nec- 
essary to look at these matters more carefully. The use of the 
perturbation technique can only be justifled retrospectively 
by comparing corrections to the uncorrelated propagation 
of the quantities involved. Here, T^ indeed turned out to be 
not more than about 0.5 of the zeroth-order quantity. If, for 
a difi^erent problem, this is not the case, one would have to 
sum the entire series 13211 or seek other solution methods. 

We have also only considered a point-like source and 
constant density of the deflectors. It does not pose any con- 
ceptual difficulties to perform the same calculations for a 
different source profile and deflectors density. However, such 
calculations are likely to require considerable numerical ef- 

Apart from these approximations, we have also made 
certain assumptions about different parameters of the lens- 
ing population and most of our numerical predictions would 
change if very different parameters were considered. How- 
ever, there are some generic features which do not suffer from 
these uncertainties, the most robust of them are the shape 
of the power spectrum broadening function ■!/'(Atj) and the 
autocorrelation function correction V'(i)- In addition, due to 
the saturation of the variation amplitude at relatively low 
q ~ 10, it is not very sensitive even to considerable changes 
in v, po or m. 

To sum up, we have shown that microlensing needs to be 
taken into account in the interpretation of the power spectra 
and variability amplitudes of various sources. These predic- 
tions are well substantiated in the case of the Galaxy, and 
therefore timing properties of pulsars and X-ray binaries has 
the potential to reveal properties of the microlensing pop- 
ulation of the Galaxy. For cosmological distances, our pre- 
dictions are weaker because the knowledge of the microlens- 
ing population is less certain and the model of uniformly 
distributed lenses is unlikely to be applicable in this case. 
When this difficulty is overcome, one can attempt to study 
the cosmological microlensing population via this method. 

Acknowledgements. AVT is supported by IPRS and 
IPA from the University of Sydney. We thank Rodrigo Gil- 
Merino for useful comments regarding the manuscript and 
the referee for suggestions which lead to a considerable clar- 
ification of the presentation. AVT also acknowledges owing 
drinks to J. Berian James for the effort he put into the work 
on the English of this paper. 


Abramowitz M., Stegun I.A., 1972, Handbook of mathe- 
matical functions, 10th edition, 503-536 

Alcock C. et al., 2000, ApJ, 542, 281 

AmanuUah R., Mortsell E., Goobar A., 2003, A&A, 397, 

Astier P. et al., 2005, A&A accepted l |astro-ph/0510447| l 

BahcaU N.A., Siang P.O., 1996, ApJ, 462, L49 

Binney J., Tremain S., 1988, Galactic Dynamics, Princeton 

University Press 
Bregman J.N., 1990, Astron. and Astrop.Rev., 2, 125 
Chandrasekhar S., 1943, Rev. Mod. Phys., 15, 1 
Deguchi S., Watson W.D., 1987, Phys. Rev. Letters, 59(24), 

de Vries W.H., Becker R.H., White R.L., Loomis C., 2005, 

AJ, 129, 615 
Edelson R., Kirpal N., 1999, ApJ, 514, 682 
Evans N.W., Belokurov V.A., 2004, preprint 

(astro-ph/04 11222) 
Fiore P., Laor A., Elvis M., Nicastro P., Giallonco E., 1998, 

ApJ, 503, 607 
Giveon U., Maoz D., Kaspi Sh., Netzer H., Smith P.S., 

1999, MNRAS, 306, 637 
Gradshteyn I.S., Ryzhik I.M., 1994, Table of integrals, sums 

and products, 5th edition 
Hafizi M., De Paolis F., Ingrosso G., Nucita A.A., 2004, 

Int.J.Mod.Phys.D 13, 1831 
Han C., Gould A., 2003, ApJ, 592 172 
Hasinger G., van der Klis M., 1989, A&A, 225, 79 
Hawkins M.R.S., 1996, MNRAS, 278, 787 
Hellier C., 2004, in Tovmassian G., Sion E.M., eds. Com- 
pact Binaries in the Galaxy and Beyond, Rev. Mex. de 

Astr., Vol. 20. 
Holberg J.B., Oswalt T.D., Sion E.M., 2002, ApJ, 571, 512 
Holz D.E., Linder E.V., 2005, ApJ, 631, 678 
Hui L., Green P.B., 2005, preprint (astro-ph/0512159) 
Katz N., Balbus S., Paczyhski B., 1986, ApJ, 306, 2 
Kayser R., Refsdal S., StabeU R., 1986, A&A, 166, 36 
Kluzniak W., Lasota J. -P., Abramowicz M. A., Warner B., 

2005, A&A, 440, L25 
Lahav O., Liddle A.R., Eidelman S., 2004, Phys..Lett. B 

592, 1 
Lasserre T. et al., 2000, A&A, 355, L39 
Lazio T.J.W., Cordes J.M. de Bruyn A.G., Macquart J.-P., 

2004, New Astron. Rec, 48, 1439 
Lewis G.F., Miralda-Escude J., Richardson D.C., Wambs- 

ganss J., 1993, MNRAS, 261, 647 
Mauche C.W., 2002, ApJ, 580, 423 
Neindorf B., 2003, A & A, 404, 83 
Oppenheimer B.R., Hambly N.C., Digby A.P., 

Hodgkin S.T., Saumon D., 2001, Science, 292, 698 
Paczyiiski B., 1986, ApJ,304, 1 
Paltani S., Courvoisier T.J.-L., Blecha A., Bratschi P., 

1997, A&A, 327, 539 
Perlmutter S. et al., 1999, ApJ, 517, 565 
Psaltis D., 2001, Adv. in Sp. Res., 28, 481 
Riess A.G. et al., 1998, AJ, 116, 1009 
Riess A.G. et al., 2004, ApJ, 607, 665 
Sadoulet B., 1999, Rev. Mod. Phys., 71, 197 
Scalo J., Elmegreen B.J., 2004, Ann. Rev. A&A, 42, 275 
Schneider P., Ehlers J., Falco E.E., 1992, Gravitational 

lenses. Springer 
Schneider P., Weiss A., 1987, A&A, 164, 237 
Seitz C., Schneider P., 1994, A&A, 288, 1 
Seitz C., Wambsganss J., Schneider P., 1994, A&A, 288, 19 
Shapiro S.L., Teukolsky S.A., 1983, Black Holes, White 

Dwarfs and Neutron Stars, Wiley 
Tonry J.T. et al. 2003, ApJ, 594, 1 

14 Tuntsov & Lewis 

Tormen G., Moscardini L., Lucchin F., Matarrese S., 1993, 
ApJ, 411, 16 

Tuntsov A.V., Lewis G.F., 2006, MNRAS, this volume (Pa- 
per I) 

van der Klis M., 2006, in Lewin W., van der Klis M., 
eds. Compact Stellar X-ray Sources, Cambridge Univer- 
sity Press 

van der Klis M., 2005, Astron. Nachr. 326, 798 

Wambsganss J., Witt H.J., Schneider P., 1992 A&A, 258, 

Warner B., 2004, PASP, 116, 115 

Warner B., Woudt P.A., 2005, in Hameury J.-M., Lasota J.- 
P. , eds. The Astrophysics of Cataclysmic Variables and 
Related Objects, ASP Conf. Series, Vol. 330. 

Witt H.J., 1993, ApJ, 402, 530 

Wood-Vasey W.M. et al., 2004, New Astr.Rev., 48, 637 

Zackrisson E., Bergvall N., Marquart T., Helbig P., 2003, 
A&A, 408, 17 


A complete eigensystem of solutions to the propagation 
problem for one-ray statistics was found in Paper 1: 

exp [— i (k • r + ujt)] 


M^Lm = c;™(t) 


Ak^nm = ^t;— k^ +i-yiK(cJ) (1 + \n\ +2m) . 
2lu c 

where {x '■ ''' ~ t^x) 

, exp[— i(nx + ^k ■ r)] 



^nm \'^ ) — 



( {/iK(a 




X exp 




The 'oscillator constant' n in the above expression is defined 

Kic) ^ :^\\c\\{f3') 


and, according to our convention, when a square or quartic 
root of in is calculated, one takes the value in the complex 
plane closest to the positive ray of the real axis. The conju- 
gate basis Af^J^^ can be obtained by complex-conjugating 
Cnm(T) function only. See Paper I for details. 

All possible combinations of the form 

Ma{ri,r2,ti,t2,Ti,T2) = M^^^ (n , ti , Ti) M^^' (ra, ^2 , Ta) 
and associated sums \a = Ajj^ +^a2 with a composite index 
a = (ki,2, ^ =F '^/2, m.i^2,ni^2) form a dense set of solutions 
to the unperturbed problem 1251 1 with zero correlation 
term. In this basis, the initial conditions for a synchronous 
isotropic source of the form I12II can be written down as: 


16nc^ g(ki)ff(k2)P(a;) 


X exp 

- (ki - k?) D. 


Q-f 1 


ijjDs\J n({jS) 



In order to obtain the expectation value of the flux second 
moment one needs to project the intensity moment vector 
71/2 onto T = 0: 

F2 = M2|^^ ^^0 


As explained in the main text, this projection leaves intact 
the spatial and temporal sectors and therefore for ki,2 = 
K =p k/2, ijj\^2 = f2 =F t^/Z, Fourier modes of the flux are 
given in Ma basis as a sum 

F2(K,k,n,c^) = ^c4D)e^"^CaU,,,.o (B2) 

with respect to angular modes only. 

In particular, since C^ranij = 0) = S^d \JiK.{uj)/-K, one 
has for fi = 0: 

F2{K,k,Lu;D) = 

71" Z — / 

mi. 2 


Note, that ni,2 7^ can be ingnored. 

These moments can be analogously split into the zeroth- 
order contribution 

F2(K,k,(^;D) = } Ck 

TV ^ — ' 



and a term due to the correlated deflections (see II36II ') 
Fi{K,k,u;-D)= (B5) 



where b = (K, k, w, 7711,2). 

The summation above is, in general, rather difficult 
to perform. However, significant simplification occurs for 
an isotropic and synchronous source of form (IA6II . First 
of all, initial conditions for isotropic sources are propor- 
tional to 5^,, hence one does not need to perform summa- 
tion in n' (these components do develop for D > but do 
not contribute to the fiux, according to JB3l l). Second, the 
set cjj in this case is completely defined by it spatial pro- 
file g(K' -'k'/2)g{K' + k'/2) and power spectrum P{uj'); 
hence, due to the linearity of the problem there should be 
coefficients Fad'{D) such that: 

F(K, \i,u;D)=J2 P^d' (^)^(^')fl(k'i)5(k': 


with d — (K,k,(jj); here we use the summation instead of 
the integration sign to avoid unnecessarily long notation. 
Noting that for sources described by IA6II 

e^'""!"^^ =e^'*^,7,(D)™i (773(1?))'"^ (B7) 


cLi^, = p('^)ff(ki)<7(k2)c«rHr)'^ 


rid{D) = exp 

-7^(1 + i)D 



a + 1 

Microlensing in phase space 11: Correlations analysis 15 








it is easy to see that in the zeroth and first order the coeffi- 
cients Fdd' are given as 



TV ^ — ' 

1 _ \/M^) 
-fdd' — Cjj/ 

J 7-i' ->>rf(D-D') A,/D' 

dD e ''^ 'e <* 


E/^ (^ l\vn-. / /-* /*\m^ 7n-\ / *\m.9 

where r] = rid{D — D') and -q' — rid'{D') and the sum is over 
w-i,2,?Tii,2 running from to infinity. 

Summation in lB13t can be done immediateiy using 

X; z"' = 1/(1 - z) vaiid "iz e Z : \z\ < 1 and gives 

F'L'iD)^ (B15) 

d 16c 

,2 exp 

-i2£ (K • k) {Ds +D)-^ V^^ 

a + 1 — (q — 1) exp 

^^V^H(l + i)^ 

In order to evaiuate sum in lB14t we recali that Co is a 
function of r and t only and therefore the matrix elements 

Cw = -^||/:||(/3') (Aft, n^ Cora Mfc) 


a'^.t- t (^ \r l^ exp(iQi-Ti) 

X / d T2C+ij/2,0,m2(T"2)C+^'/2,0,m' (t2) 

exp(iQ2 ■ T2) 

xr/ Co(k - k',a; - uj')t2 ( --||/:| |(/3")5(K - K'] 




K =F k/2 K' =F k72 


and C'o(Ak, Alj) is defined by 12111 . Functions ^u,o,m are 
given by (IA4II : 

Ca;,0,m(T) = CS-'^m f •\/iK(c<j)7 


with the coefficient CZ = y 2^iK(Lj) exp — |y'iK(a7)r^ 

dependent of m. 

Therefore, one can change the order of summation with 
respect to m and integration with respect to r in 1B14I I and 
make use of the Laguerre polynomials generating function: 







to obtain for the sum in the second line of IIB141 

J(D') ^-'iu\\C\\{lf)^K(Lu)K{Lo')5{K-K!) (B20) 

exp [-7riV2 - -y'T^^] 

j2 j2 
d Tid T2 


X (27r)-^c'«-^i'n^Co (k - k', a; - J) rae'^'^-^^^ 


^ , /IT-TTTi+i!/ 

7 = (1 - i) vM^Y^ + vM^j 



If one now denotes yii,2 the directions of Qi,2 
<y5i,2 : Qi,2 = Qi,2e^i2, angular integration in l)B20^ 
yields a product of two first-order Bessel functions and 
radial integration in r can be performed using iden- 
tity fp°°drTVi(QT)exp(-7 rV2) = g/72 exp(-gV27) 
JGradshtevn fc Rvzhik^ll994^ . Thus we arrive at 

I{D') = AvyfK,{uj)K(Lu')5 (K - K') C (Ak, Alj) (B22) 

Q1Q2 exp [— Q?/27 — Q2/27*] cos(v3i — a) cos(q — 1^2) 

72(1 - r,)(l - ^V) 7*2(1 _ ^.)(1 _ eV*) 


, _ ,^<:Jr 16c^ exp[-ig£(K'.kOO.] 
-f'dd'(W) - \/ -7-77 -777^ n , ,.,12 (B23) 

fi:(a;') Lj'2i32 

1 + a- 

I AD' I{D') exp [Xd{D - D') + Ad-D'] 


In this appendix we calculate quantities 13911 and the line 
profile function %j) in the simple case of a point-like source at 
the origin of the source plane (?(k) = 1 as seen by a single 
observer at the origin of its plane ri,2 = 0. 

Using IB15II one obtains, in this situation (care must be 
taken when calculating the limit of k — > because a tends 
to infinity in this limit): 





Df {Ds + DY 

a + 1 — (q — 1) exp 

-fy^(l + i)Z) 


{Ds + D) 


Therefore, in the leading order, the line profile function 1/) is 
given by the expression 1)41^ . 

The last line of lUH also shows that the zeroth-order 
correction to the free propagation of fiux is of the second 
order in k (and therefore v); as shown below, correction in 
the first order of perturbation due to correlated defiections 
is O {y/K\ and therefore represents the dominant term. 

To perform integration 14H for the first-order correc- 
tion, we first untilize the 5-function in K — K' from 1B23I I 
and then, taking into account that C depends on Ak = k— k' 
only, perform a change of integration variables (K, k, k') ^ 

16 Tuntsov & Lewis 

(Qi, Q2, Ak) according to 



k = 




-^(Qi+Q.) + ^Ak 



The Jacobian determinant of this transformation 

£'(K,k,k') _ / uuj' ^ " 
D(Qi,Q2,Ak) ^ \^2cAcj 

is independent of integration variables and thus appears as 
a factor in front of the integrals, and since this is a linear 
transformation Qi,2,Ak independently run through entire 
VJ^ planes. 

Then, we may rewrite the first-order correction 14H as 

i:^{LO,J ,D) = 



{Ds + DY 

■ exp 


4i^A/K(a;)K(a)')exp -D' {^^J k{uj') - ^a/k(cj)) 



d^AkC(Ak,Aa;) / dQi dQ 




27 27* 

X / d(y9l d(y92 


COs((^l — Of) COs(</32 — Ct) 

1— (K ■ k) (D - 7:)') 


X exp -1— {iS^- k.}[U - U ) -\— (K' ■ k') (Ds + D') 

here we deliberately postpone some obvious cancellations to 
keep the asymptotic behaviour clear and provide for lens 
population parameters dependence on D. 

Using IIC2II we rewrite the argument of the exponent in 
the last line in terms of Qi,2 and Ak, angular integration 
yields a product of two Bessel functions and the identity 
used in deriving lB22t brings the integral with respect to 
Qi,2 to the following form: 

' p^Afc^ 

u(Afc) = 


\^ uju' ) (Ds+D)- 
with (|7|^ = Im^7 + Re^7) 

■ exp \—p Re7 Afc 1 

p={Ds + Df I 

2Alj^ ^ ^ 
h 5Im7 + 5 Re 7 




Function ii(Afc) is independent of its argument vector 
direction, and integration of (^(Ak, Ao;) with respect to the 
angular coordinate produces Infeld (modified Bessel) func- 
tion 27r/o(uALi;/cr^Afc). Therefore, one has for the integral 
with respect to Ak: 


dAfc exp 




2 \ cr2Afc2 

2v2Tvm p e 2a 
a (Ds + Df 

+ p{j + 7*)AA: 


-'"(^)h(T-') -'•(*?-' 

The exponential factor decreases rapidly at Afc tending to 
both zero and infinity thus providing necessary convergence. 

It is reasonable now to introduce some approximations. 
It was s hown in Paper I that in physically reasonable situa- 
tions ^yK{Lo) <C 1. We can use this fact and expand 77, ^, r;' 
etc. in its power series retaining only the few leading terms 
as necessary. 

First of all, we have for the denominator in I{D') 

f ^2K{Lj){^^{D-D')y +o{k) 


\l-iri'\^=2K{Lu')(^{Ds + D'] 

and similarly 

1 lu'''DU{lu') 

l + a|2 ^ 2c2 





For 7 the first term in expansion is purely geometrical and 
imaginary while the first real term is proportional to «:: 



ujuj'{D- D'){Ds + D') 



Ds + D 


{Ds + D){Ds + D') 


+ 0(«'), 

where we have used u)k(u}) = u)'k,{u)') from IIA5II . 

Expanding Auj^ and using the above formulae one finds 
that at this level of approximation p is purely geometrical 
in the leading order: 

p = 

(D, + i?)2 + 0(^2 


with dimensionless distance coordinate t = {Ds + D')/{Ds + 
D) introduced; as D' changes from Ds to D, the correspond- 
ing range for t is from s = Ds/{Ds -I- D) to 1. 

Now, assuming that parameters v, m, a, v are indepen- 
dent of D', we use HASH and {/3^) = 27r||/:||"^m^ In A, intro- 
duce the frequency difference scale 




dimensionless variables 
x = ^Afc, 




and an amplitude parameter 

2ljk{uj) 2 (Ds + Df'^ 






m /87rlnAi/(Ds + -D) 

we write down the following approximation for the ra- 
tio llC4|l : 


12 exp [-uV2o-^] 

2tt In A Aluo 


Microlensing in phase space 11: Correlations analysis 17 

Figure CI. Auxiliary function w{s,£) defined by the inte- 
gral iC20i . Behaviour of «) as a function of J is displayed for 
five different values of s - 0, 0.2, 0.4, 0.6 and 0.8; w(s,i) oc ^"^ 
for large 1^. 

Xq I <it <ixt (1 — t) exp 

2 2 

-— -(1-i) {t -s ) 

X exp 



^0 f-^) [Jo (x) - Jo {Ax)f + o{V^). 

- as expected, it is 0{y^). 

From this point, integration needs to be numerical, and 
is, perhaps, most easily done by first calculating the integral 
of the function on the second line of IClSt with respect to t 
and then performing calculation in x. This approach yields 
the expression 14411 with the scaled line profile function 




q I dxw {s,xq) exp 

defined through an auxiliary function 

w{s,^)= dtt^{l~tfexp 


Jof-^) [Jo {x) - Jo {Ax)f 
\a X J 



The behaviour of w(s, ^) as a function of ^ is shown in Fig- 
ure lUll for a few values of s; w(s,5) oc ^~^ for ^ S> 1. Since 
A S> 1 we can safely neglect the second Bessel function in 
the square brackets when evaluating (IC19II .